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The long-standing problem of finding a closed formula for the relative entropy of entanglement 
(REE) for two qubits is addressed. A compact-form solution to the inverse problem, which charac- 
terizes an entangled state for a given closest separable state, is obtained. Analysis of the formula for 
a large class of entangled states strongly suggests that a compact analytical solution of the original 
problem, which corresponds to finding the closest separable state for a given entangled state, can be 
given only in some special cases. A few applications of the compact-form formula are given to show 
additivity of the REE, to relate the REE with the Rains upper bound for distillable entanglement, 
and to show that a Bell state does not have a unique closest separable state. 
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I. INTRODUCTION 

The relative entropy of entanglement (REE) is as an 
entanglement measure quantifying how much a given en- 
tangled state can be distinguished operationally from the 
set of separable states or those with positive partial trans- 
position (PPT) Q: 



Er(p) = min^ eV S(p\\a') = S(p\\a), 



where V denotes a set of separable states or PPT states, 
S is a quasidistance measure usually chosen to be the 
quantum relative entropy, 5(/9||cr') = tr (p\gp — p\ga'), 
an analog of the classical Kullback-Leibler divergence. A 
state a on the boundary of separable states is called the 
closest separable state (CSS) or the closest PPT state. 
Various properties of the REE have already been de- 
scribed (see, e.g., 0, S ! I, 0, [EM 1 I, GSM GJ, Dl 

[2 EE EH and for a review se e 1171 ). But there is still 
an open fundamental problem [181 ]: the task of finding a 
closed explicit formula for the REE for a given two-qubit 
state corresponding to a solution of the convex optimiza- 
tion problem for the REE, or, briefly, of finding the CSS 
a for a given entangled state p. 

Here, we give a compact-form solution of the closely 
related problem of finding entangled states p and their 
REE for a given CSS a. Our formula is derived from 
the results obtained by one of us in Ref. EH]. We also 
demonstrate the intrinsic difficulty in solving the original 
problem. In addition, we apply our formula to relate 
the REE with the Rains upper bound for the distillable 
entanglement, and to show additivity of the REE and 
nonuniqueness of the CSS for a Bell state. 



II. A CLOSED FORMULA FOR THE REE 

Let us consider an entangled two-qubit state p with its 
CSS er, and assume that a is full rank. Note that the case 
includes the arbitrary full-rank entangled state p, since 
rank(cr) > rank(p) must hold so that S(p\ \a) is finite, and 



hence the CSS a for full-rank p is always full rank. Since 
a is an edge state then its partial transposition a r is rank 
deficient, i.e., rank (cr r ) = 3. Let |0) be the kernel (or 
null space) of <7 r , i.e., an eigenstate of o~ r corresponding 
to zero eigenvalue, 



= 0. 



(2) 



Moreover, let \i) and Xi be eigenstates and eigenvalues of 
(1) a, respectively. Using the kernel \(f>), the formula in Ref. 



[131 ] can be rewritten in the following simple form: 

xG(a), (3) 



where 



G(a) 



G, 



E 



A, 

IgAi-lg \j 



for Xi = Xj , 
for A l ^ Xj , 



(4) 



(5) 



and x > 0. All p obtained from Eq. ([3]) for x max > x > 
have a as their CSS, where x max is the threshold for 
p > 0. This is a unique solution of all extremal conditions 
in two qubits as shown later. The same relation holds 
for entangled states and the closest PPT state even in 
higher-dimensional systems as long as the (PPT) entan- 
glement witness (EW) Z such that tr Za = is uniquely 
determined as (cj)\) r . Moreover, from Eqs. ©-((SI) we 
have 



(i|p|i) = A i [l- a; <i|(|0)^|) r |i)]. 
Therefore, the REE can be rewritten as 
E R (p) = trplgp - trplgcr 

= tTplgp-y~){i\p\i)lg\i 



(6) 



= tr p lg p — tr a lg a 

+x^(z|(|0)(0|) r | i )A l IgXi 

i 

= S(a)~S(p) + xtT[(\cj } ){cb\) r a\ga], 



(7) 
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where S(-) is the von Neumann entropy. In any case, 
however, Eq. © should be conversely solved with respect 
to a to obtain the true closed formula for the REE as a 
solution to Eisert's problem [l8|. In the following we 
show that the inversion is possible in some special cases 
but rather not in general. 

Our formula can be derived from the results obtained 
in Ref. [HI but the derivation is rather lengthy. For- 
tunately, we can also derive it in a compact and elegant 
way using a different approach based on the result shown 
in Ref. that the operator of 



Z 



1 



P- 



a + z a + z 



-dz 



(8) 



must be an entan glem ent witness. Let us briefly repeat 
the proof in Ref. [1 91 ] for our later convenience. Since a 
is the CSS for p, the inequality 



S(P\\(1 



e a 



S(p\\v)>0 



(9) 



must hold for every separable state a' and < e < 1. 
Using the expansion of 



\g(X + eY) = lg(X) + e f 
Jo 



1 1 



Y——dz + 0(e 2 ), 



X + z X + z 



we have 



S(p\\(l~e)a + ea')-S(p\\a) 

f°° 1 1 

= etrp / O-o-') dz + 0(e 2 ) 

J a + z a + z 

r°° i i 

= e[trp-trer' / p dz] + 0(e 2 ) 

J Q a + z a + z 

= etr Za' + 0(e 2 ), (10) 

and hence etr Za' + 0(e 2 ) > must hold for arbitrary 
small e > if a' is a separable state. This implies that 
tr Za 1 > must hold for every separable state a', and 
therefore Z must be an EW (l9| . 

Now suppose that a and p are two-qubit states. It was 
shown in Ref. [2(| that an EW in two qubits must be 
decomposable, since there are no PPT entangled states 
in two qubits. Therefore, Z must be a decomposable 
EW, and hence Z = P + Q T , where P and Q are positive 
operators [2(j. Moreover, 



a + z a + z 



: v Za = 1 r a- - /> / n d : = tr [a — p) = 0, 

(11) 

and as a result tr Za = tr Pa + tr Qa L = must hold. 
Since a is full rank, and a r is positive and rank 3, the 
solution is uniquely determined (leaving out the normal- 
ization of x > 0) as P — and Q = x\(j>){<fi\. As a result, 
x(\4>) (4>\) T = Z holds for full rank a. The integral in Z 
can be performed using the eigenstates \i) for a such that 



*<i|(|<^|) r |i> - 



i 



o A, + z 



'Mj) 



\j + z 



-dz 



(12) 



and finally we have Eq. ([3]). In this way, the satisfaction 
of all extremal conditions in the optimization problem for 
the REE is automatically ensured by the condition that 
Z is an EW such that tr Za = 0. Note that \cf>} is always 
entangled for full rank a [2l| (if \<j>) is not entangled, the 
full rank a is not a CSS for any entangled state). 



III. THE RAINS BOUND AND THE REE 

An upper bound for distillablc entanglement intro- 
duced by Rains [6(] is defined as 



R(p) = miE L[S(p||r')+lgtr|r' r |], 



(13) 



where minimization is taken over all states including en- 
tangled states, and hence R(p) < Er(p) follows from the 
definition. Here, let us apply the technique as used in 
the previous section to the optimization problem for the 
Rains bound. It was shown in Ref. [ll| that the problem 
is reduced to 



R(p) = mm S(p\\T>), 

t'>0 



(14) 



where the minimization is taken over unnormalized states 
subject to tr|r' r | < 1. This is a convex optimization 
problem because tr It^I < 1 for To — pr\ + (1 — p)T2 and 
tr |rf |, tr \t% \ < 1 [D]]. Suppose that p is full rank and r 
is an optimal unnormalized state. Hence, r is full rank, 



tr r 1 



1, and 



5(p||(l-e)r + er')-5(p||r)>Q (15) 

for every r' such that tr |r' r | < 1 and < e < 1. Let r' be 
a normalized separable state, i.e., r' r > and tr |r' r | = 
tr t' = 1 . Using the expansion of the logarithmic function 
for e — > +0 as in the previous section, we have 



S(p\\(l-e)r 

- etr \l - 

Jo T 
= etr Z r t' > 0, 



S(p\\r) 

— p dz] r . 



(16) 



and hence Zr must be again an EW. Note however that 
tr Zrt = tr t — tr p < in this case, contrary to Eq. (|ll|l . 
because r is unnormalized so that tr |r r | = 1. 

Let us then consider the case where p is a two-qubit 
state, and suppose that the optimal two-qubit state r is 
entangled. Since the partial transposition of a two-qubit 
state has only one negative eigenvalue [9j] , T r is expressed 
such that 



= (i-M)n 



(17) 



where n > 0, IL\ip) = 0, and p > 0. Moreover, trll = 1 
so that tr |r r | = 1. For a small deviation of p — > (l + 5)p, 
we have 

S(p\\r) - S(p\\r)-5p[trZ R (U r + (\i;)(i ! \) r )-2] 
= S( P \\t) - S [tr Z R (I? -t)-2u] 
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S( P \\t) - Six Z R I? , 



(18) 



where tr Zrt = trr — 1 = — 2/z was used. Since r 
must satisfy the extremal condition with respect to p, 
tr ZuII r = must hold. Moreover, since x is a two- 
qubit entangled state, II is rank-3 and II r is positive 
definite and as a result the EW Zr is uniquely de- 
termined (including the normalization in this case) as 
Zr = 2(\tp) (ip\) r , where tr Zrt = — 2p and \x ^ were 
taken into account. This implies that 



i 



T + Z T + Z 



dz, 



(19) 



but this cannot be satisfied because the right-hand side 
is positive definite for full-rank p and r, while the left- 
hand side is not for any \ip). Therefore, the optimal state 
r must not be entangled, and the optimization in R{p) 
is achieved by a separable state. The same discussion 
also holds for low-rank p, because R{p) is a continuous 
function [ll|. It is then concluded that R(p) = Er{p) 
for every two-qubit state, and our compact-form formula 
also holds for the Rains bound. 

Note that the Rains bound is strictly smaller than the 
REE for the Werner state in higher-dimensional systems 
[H> @> but such a disagreement does not occur in two 
qubits as shown above. 



IV. ADDITIVITY OF THE REE 



An asymptotic REE defined as 



lim -E R (p® n ) 

n—*oo ft 



(20) 



satisfies E^(p) < Er(p) from the definition. The equal- 
ity holds if Er(p) is weakly additive, but this is not the 
case in general Here, let us briefly investigate the 
additivity using our compact-form formula. 

In Ref. [26], it was shown that Er(p) for p such that 
[p, a] — is weakly additive if (per^ 1 ) 1 " > — 1. From Eq. 

m, 



(i\(lg Xip- pig Xj)\j) 

= -^|[A 4 (|0)(0|) r 



(1^1)%] \j) (21) 



must hold for all i and j, and hence we have [lgc, p] = 
—x[a, (\cf))(<f)\f]. This implies that [p, a] = if and only 
if [a, (|0)(^|) r ] = 0. Therefore, (\4>}(4>\) r must be diago- 
nalizcd in terms of the eigenstates of a so that [p, er] = 0, 
and the compact-form formula in this case is much sim- 
plified as 



p = a- X (\4>)(<P\f<T. 



(22) 



Let po > 1/2 be the maximal Schmidt coefficient of \4>). 
Since the largest eigenvalue of (\<j)} (4>\) T is po, the range 
of x must satisfy < x < l/po < 2, so that p > 0. 
As a result, (pcr _1 ) r = I — x\<fr)((j>\ > — 1 always holds. 
Therefore, it is concluded that Er(p) is weakly additive 



and E^(p) — Er(p) for every two-qubit state such that 
[p,a]=0. 

Moreover, it was shown in Ref. [26| that Er(p) for p 
such that [p, a] — is strongly additive, namely Er(p <g> 
t) = Er(p) + Er(t) for an arbitrary r, if (pcr _1 ) r > 0. 
In the same way as above, it is found that the condition is 
satisfied if x < 1, and therefore Er(p) is strongly additive 
for every two-qubit state such that [p, a] = and x < 1 . 



V. FORMULA APPLICATIONS FOR 
FULL-RANK CSS 

All the examples of arbitrary rank states p with their 
CSS a found by us in the literature [E 0, i, i, M, OH EE 
[l6j , can easily be explained using our formula. The pro- 
cedure can be summarized as follows: choose a full-rank 
matrix a, calculate its partial transposition to get a T , 
find a condition on its elements for which er r is rank de- 
ficient (and so becomes a CSS), calculate p = a — xG(a), 
if required take a limit of some elements to diminish the 
rank of a and p, and finally find an inverse relation to 
express the elements of a in terms of those of p. 

For example, let us analyze a full-rank state a = 
J2i=i diagonal in the Bell basis {|/3i)}- The 

eigenvalues of a r are A, = \ — Ri. Thus, e.g., by setting 
Ai = 0, a becomes a CSS. By noting that the kernel \4>) 
is a Bell state and applying it to Eq. one gets a 

Bell-diagonal entangled state pbd = Yli r i\Pi)(fli\> where 
n = (2 + x)/4 and otherwise ^ = Ri(l— x/2). By invert- 
ing the latter equation, one gets the well-known formula 

Q 



n) 



Y,rmm\- (23) 



This is the CSS of an arbitrary Bell-diagonal state pbd = 
Ei=i r i\Pi){Pi\ assuming n > 1/2. 

As another example which, to our knowledge, has not 
been discussed in the literature, let us analyze a two- 
qubit state of the following form: 



ARi 

i? 2 Y 

Y i? 3 

V Rij 



(24) 



This state is the CSS if its partial transposition <j\ is rank 
3, which implies that Y = \JR\R±e lip (in the following 
we set ip — 0), while the requirement of positivity of 
the density operator ctz implies that R2R3 > Y 2 . Thus, 
Pz satisfies the condition R2R3 > R1R4, and (f2"4")l can 
compactly be given by 

a z = AA(|i/')(V'l) r +i?2|01)(01| + i?3|10)(10|, (25) 

where \i/>) = Af- 1/2 (y/R[\00) + VR~l\ll)) and N = Ri + 
i?4. The eigenvalues of e*"z are {Xi)i — (Ri, i?4, A + , A_), 
where A± = 5(^2 + R3 ± z) with the auxiliary func- 
tion z = J (R 2 — R3) 2 + 4Y 2 . The corresponding eigen- 
vectors are |Ai) = |00), |A 2 ) = |11), and |A±) = 



4 



A/±[(A± - i? 3 )|01) + Y\W)] with normalizations A/± = 
[(A± - R3) 2 + Y 2 ]- 1 / 2 . One finds the kernel \<j>) of a\ 
and then 

(|(/»)(0|) r = (i?.! + J R 4 )- 1 (i? 4 |oo)(oo| - r|oi)(io| 

-y|10)(01| + i? 1 |ll)(ll|). (26) 
Thus, according to ([4]), we find that 



G(az) 



where 



Rl = i?4 



(Ri \ 

i? 2 Y 

y i? 3 

\ i? 4 / 



1- 



(27) 



R 2 

R 3 
Y 



Ri + i?4 

-2F 2 rf[(i? 2 - i? 3 )(z + R 2 L) + 2Y 2 L], 
—2R\ — i?2, 

Yd[2Y 2 (R 2 + R 3 )L- (R 2 -R 3 ) 2 z] (28) 



are given in terms of 



z = ^{R 2 -R 3 ) 2 + 4i?ii?4, 

L = ln(E 2 + i? 3 - 2) - m(fl 2 + #3 + z), (29) 

and l/d= (Ri + R 4 )z 2 L. 

Thus, according to ©, the entangled states 



Pz = <?z - xG(a z ) 



(30) 



have the same CSS rrz- This is an important example in 
our analysis as all low-rank states discussed in Sec. VI are 
special cases of (13T7|) . Note that it is required to assume 
x < ^max = {R1+R4)/ Ri to ensure that (R4 — XR4) > 0. 
Assuming for simplicity that Ri > R4, the analogous 
conditions for (Ri — xRi) > with i — 1,2,3 are also 
satisfied for x < x' mliK . On the other hand, the condition 
(R 2 ~xR 2 )(R 3 — xR 3 ) > (Y — xY) 2 , which is also implied 
by the positivity of pz, restricts x to be smaller than 
< ax = / - v// 2 -4AA/(2A), where A = R 2 R 3 - Y 2 , 
A = R2R3 - Y 2 , and / = R 2 R 3 + R 2 R 3 - 2YY. Thus, 
(|30"f is defined for < x < x max = minjx^^, x^ ax }. The 
problem of expressing az in terms of pz will be addressed 
in the following section. 



VI. FORMULA APPLICATIONS FOR 
LOWER-RANK CSS 

Our compact-form formula can also be applied for 
lower-rank CSSs a in two approaches: directly for some 
special states and indirectly for arbitrary states. 

To justify a direct application, the following conditions 
should be satisfied: (i) There must exist a full-rank edge 
state cr' in the vicinity of a. If it is certain that a is a CSS 
for some p, this condition is trivial. However, if we do not 
know whether or not a can be a CSS for some p, this is 



not trivial, (ii) Let be the kernel of cr' r , i.e., a lT \4>') — 
0. Then, \<fi') must be entangled. However, when \4>'} 
is not entangled, we merely cannot find any entangled 
p' > for a' by the compact-form formula, and hence this 
condition is not so important, (iii) There must exist a 
sequence such that a' — > a and |</>') — > \<j>) simultaneously. 
This condition seems to severely constrain the choice of 
\(j)) in the case of rank(c r ) = 2. 

Thus, our formula can be applied directly to the rank-2 
Horodecki state defined for p £ (0, 1) by [22| 



ph = p|^ (±) )(V' (±) I + (i-p)|oo)(oo|, 



(±) 



(31) 



which is a mixture of a Bell state IV'^) = (|01) ± 
\10})/^/2 and a separable state orthogonal to it. It is 
worth noting that the Horodecki state is extremal in the 
sense that it minimizes the REE for a given concurrence 
Q, negativity (i.e., a measure of the PPT entanglement 
cost) for a given concurrence @, fidelity (i.e., maximal 
singlet fraction) for a given concurrence (> 1/3), and 
negativity [> (\fh — 2)/3] [2||. The state also satisfies 
some extremal conditions for the REE for a given neg- 
ativity 45]. The CSS for p H derived from ©-© reads 
as 



oh 



g' 2 |00)(00| + 2 P y|V (±) )(V (±) l +P' 2 |11)(H|, (32) 



where p' — p/2 and q' = 1 — p', in agreement with the 
known result derived in another way [2|. Note that, al- 
though rank (o^) = 3, cth is not full rank. Fortunately, 
the above-mentioned conditions necessary for the direct 
application of the compact-form formula are satisfied. 

In the second more general approach, one can apply 
our formula for arbitrary lower-rank states in a limiting 
sequence from a full-rank state because the REE is a 
continuous function. For example of such application of 
our formula for lower-rank states we will analyze pure 
states and the rank-2 Vedral-Plenio state defined by 0: 

pvp =p|V (+) )(^ (+3 | + (1 -P)|01)(01| (33) 
for < p < 1 with the corresponding CSS 



ovp 



1 - |) ]01><01| + |jl0><10], (34) 



for which rank (cvp) — rank (fTy P ) = 2. By contrast with 
the Horodecki state, ([3"3"| is a mixture of a Bell state and 
a separable state not orthogonal to it. 

To derive CSSs up and cryp, and thus to show the 
usefulness of our formula also for lower-rank states, let 
us apply state <7z, given by (|24|) . in the limiting cases. 
Namely, by assuming in (f3"0")) that x — x( nax < x"^ ax and 
i?i > i?4, one gets the following extremal state: 



Pz = Pz{x = x' max ) 



(r x 0\ 

r 2 y 

y r 3 

Vo 0/ 



(35) 
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where 



T2 = P 2 + ^(P 2 2 -P 2 P 3 + 2V 2 ) + ^i(P 2 -P 3 ), 
r 3 = 1 - n - r 2 , (36) 
2/ = ^P^i + ^)P 4 - (r 2 - R 2 )(R 2 - R 3 )\. 
In the special case of a' z for R± = R± — > 0, one gets 

<4 = (1-P 3 )|01)(01|+P 3 |10)(10|, (37) 

then = 4' - ieF(|01)(10| + |10)(01|) with Y = -(1 - 
2P 3 )/[4atanh(l — 2P3)], where atanh is the inverse hy- 
perbolic tangent. This state for x — R3/Y assuming 
R3 < 1/2 corresponds to a generalized Vedral-Plcnio 
state: 



Pgvp 



phM<tM + (i-p)|oi)<oi|, 



(38) 



where in comparison with the standard Vedral-Plenio 
state pvp, given by ([55)1 , a Bell state is replaced by a 
pure state 



|Vp) = VP|01) + e^VTTPllO), 



(39) 



for any < P < 1. For simplicity we set <p = 0. Thus, the 
CSS o-gvp is just given by Eq. ([37]) for_P 3 = p(l - P). 
By assuming P3 = p/2 and x = P3/P, one gets pyp 
and its CSS ovp, given by (|34p . in agreement with the 
solution obtained in Ref. Q. Moreover, any two-qubit 
pure state |^) = c |00) + ci|01) + c 2 |10) + c 3 |ll) can be 
transformed by local rotations into the state (f3"9"| . Thus, 
the CSS a z , given by (j37[) , also describes the CSS op for 
an arbitrary pure state \if)) as expecte d Hf . 
A mixed state introduced by Gisin [23| . 

p G = g|QO)(00|+p|Vp)<Vj»|+«|ll)<H|, (40) 

where IV'p) is given by (I39)) , < p < 1, and g = (1— p)/2. 
By contrast to the generalized Vedral-Plenio state and 
the generalized Horodecki state, PO)) is a mixture of an 
entangled pure state and two separable states orthogonal 
to it . The Gisin state is also a special case of pz assuming 
Pi — xRi = P4 — 2P4 = q. Thus, its CSS is equal to ctz, 
given by (f2~4")l for Pi = P4. It is worth noting that PD| is 
one of the simplest examples of an entangled state, which 
does not violate any Bell-type inequality (for some range 
of p for a given P) [24| . 

Assuming y 2 — r 2 r 3 , the state p' z reduces to a rank- 
2 state, which we refer to as the generalized Horodecki 
state defined as 



Pgh 



PhM(iM + (i-p)|oo)<oo|, (4i) 



where p = 1 — n and \ipp) is given by (|39[) for P = 
r 2/(l — Ti). Note that, for P = 1/2, which corresponds 
to r 2 = r 3 , pgh reduces to the standard Horodecki state 
Ph given in terms of a Bell state \ipp) — |?/^ + ^)- On the 
other hand, the state pgh for t\ = (1 + C)/2 and r 2 — 



d+ reduces to the Verstraete-Verschelde state defined for 
C < 1/3 by H| 



(1ML x 



d + f 



c 

2 



(42) 



y 000/ 



where d± = \(l-C ± VI - 2C - 3C 2 ) and C = C(p v ) 
is the concurrence, which can be expressed in terms of 
negativity N = N(p v ), as C = i[7V + ^V(4 + 5A)] 



holds. It is worth noting that state (|42ll minimizes fidelity 
for a given concurrence and negativity [25j . 

It is seen by analyzing p' z as a function of elements 
of a' z that it seems impossible to invert the general Eqs. 
(I35[) and (|36p. However, the equations can be inverted in 
some special cases. For example, by assuming R 2 = R 3 
for the state p' z , which implies r 2 — r 3 = (1 — r±)/2, one 
finds that the general relation for r 2 in Eqs. ()36|1 reduces 
to r 2 — R 2 + P4 = 1 — Pi — R 2 . Under this assumption, 
the set of equations (f36|) can be solved for {Ri} in terms 
of {ri} and y < r 2 as follows: 



P4 
R 2 

Pi 



Ar iy 2 



(1 + ri) 2 - 4y 2 ' 
P 3 = r 2 - P4, 
ri +P 4 . 



(43) 



This state, in the special case of r 2 = y, reduces to the 
standard Horodecki state pn, given by (|31[) . However, it 
does not reduce to the generalized Horodecki state pgh 
if P 7^ i, for which we can give only a formal solution 
for the REE: 

4 GH) = -H 2 (n) - ri log Pi - f_ log A_ - fl log A+, 

where H 2 (-) is the binary entropy and f± — jV±[(A± — 

P3)V^2 + ^v^]' wnne ^± an( i Mt are defined below 
Eq. (f2"5)l . A compact-form explicit formula for the REE 
for the states a' z with elements given by (|4"3"|) as well as 
for the Horodecki state or the Vedral-Plenio state is thus 
obtained. However, it seems impossible to invert (|36[) in 
order to express all {Ri} in terms of {r^} and y, even for 
the generalized Horodecki state with P, p ^ 0, ^, 1. 



VII. NONUNIQUENESS OF THE CSS FOR 
BELL STATES 



Analysis of our formula and the above examples en- 
ables us to find that for a given entangled state there 
is not always a unique CSS due to a limiting procedure. 
For this purpose let us derive CSSs for a Bell state \^ +s> ) 
from oz, given by ([24]) . First, assume that Pi = P4 = e 
and R 2 = P3 = 1/2 — e for small e > 0. By taking the 
limit of e — > 0, one gets the following CSS: 



lim a 7. 



1 



(|01)(01| + |10)(10|), (45) 
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for which G(a' Bcll ) = -±(|01)(10| + |10)(01|). By not- 
ing that x max = 2, one finds that a' BcU - x max G(<TB ell ) 
corresponds to the Bell state \ip^). The same CSS is 
obtained from special cases of the CSS for pure state and 
the Vedral-Plenio state, given by (|34| forp= 1. On the 
other hand, let us analyze crz assuming Ri — 1/4 — e for 
i = 1, ...,4. Then, in the limit of e — > 0, uz reduces to 
the CSS:' 

< u = limcrz = I (|00)<00| + 2|^ (+) )<^ (+) | + |11)<H|) , 

(46) 

for which state o BcVl — x maK G(a Ben ) also corresponds to 
the Bell state \^ + ^). Note that (|4*6")l is a special case of 
the CSS for the Horodecki state, given by (j32j) for p = 1, 
and for the Gisin state, given by (|4TJ)) for p = P = i. 

In a more general approach, let us analyze a Bcll- 
diagonal state of the form 

4 

Pbd = (1 - ke)\ ft) (ft | + A,] ft) (ft | , (47) 

i=2 

where < fcj < oo and k = k 2 + fc 3 + fc 4 such that fee < 1. 
The state in the limit of e — > reduces to the Bell state 
|ft), for which the CSS depends on {ki} as follows: 

1 1 4 

i=2 

according to (f23|) . Thus, an arbitrary Bell-diagonal state 
with one of its eigenvalues equal to 1/2 is the CSS for a 
Bell state. In special cases, OBeii goes to (|4"5j) for k 2 = 
l,fc 3 = A 4 = and to (j4l| for k 2 = 0, fc 3 = k 4 = 1 
assuming |ft, 2 ) = |V (±) > and |ft, 4 ) = |0 (±) ) = ^(|00)± 
1 11)). Other CSSs of Bell states, which are not diagonal 
in the Bell basis, can be obtained by rotating ctbcII- For 
example, the CSS 

<ell = \ (2|^ (+) )^ (+) | + MM + I^X^l) , 

(49) 

where — -^(l^ ') + |<M + ^)), is obtained by rotating 

°Beil g iven by ©. 

It is worth noting that the generated Bell-diagonal 
state p — CBcii — xG(aBcii) is independent of the pa- 
rameters {ki} in (TBoii only for x = x max , but it de- 
pends on the choice of {ki} for x < x max , although the 
largest eigenvalue of p, X\ = max(eigp), is {ki} inde- 
pendent. Thus, the REE is also independent of {ki}, as 



Er(p) = 1 - H 2 (Xi) for Ai > | and E K (p) = other- 
wise. Similarly, other entanglement measures, including 
the negativity and concurrence are {ki} independent, as 
N(p) = C(p) = 2Ai — 1. By contrast, violation of Bcll- 
inequality by p depends on all eigenvalues {A^}, so it 
depends on the choice of {ki} for ciBeii- This can be 
seen explicitly, by analyzing the Horodecki parameter M 
[271 . |28| | describing a degree of the violation of the Bell 
inequality due to Clauser, Home, Shimony, and Holt. 
VIII. CONCLUSION 

We have addressed the long-standing problem of find- 
ing a qubit formula for the relative entropy of entangle- 
ment 18] or, equivalently, of finding the CSS a for a 
given entangled state p. We have obtained a solution to 
the inverse problem by finding a compact expression for 
an entangled state p for a given CSS <r, which is a crucial 
simplification of the former solution [l3j]. The useful- 
ness of our formula can be demonstrated by finding the 
REE for some special states but also by analyzing general 
properties of the REE. Thus, we have studied (i) weak 
and strong additivity of the REE, (ii) how the REE is 
related to the Rains upper bound for the entanglement 
of distillation, and (iii) nonuniqueness of the closest sep- 
arable states for Bell states. 

All the examples of entangled states p with ana- 
lytical expression for the CSSs cr, discussed in Refs. 
0, 12, !, S, E3, EH, d, Ell , can easily be explained by us- 
ing our formula as follows. By starting from some special 
a one should generate p and then try to find an inverse 
analytical relation to express a in terms of p. Thus, for 
example, we have derived the well-known formulas for 
pure, Bell-diagonal states, the Horodecki states [H)], the 
Vedral-Plenio states 0] , and the Gisin states [23[ but also 
obtained new formulas for some other entangled states. 

We have analyzed more general states pz with elements 
{r^, which can be generated from the CSS cz with ele- 
ments {Ri} via Eqs. (|36p . The point is that, apart from 
some special cases including R 2 = R3 and R% = R4 = 0, 
the set of Eqs. (j56"| for {t^} as a function of {Ri} seem- 
ingly cannot be solved for {Ri} due to the presence of 
logarithmic functions (|29|) . Thus, we cannot express az 
in terms of pz in general. Although this is not a proof of 
impossibility, our analysis of the formula for pz in terms 
of crz strongly suggests that the inverse problem, which 
corresponds to finding a compact-form relation for a in 
terms of a given p, can be solved in some special cases 
only. 
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